arXiv:1501.05662v3 [cond-mat.dis-nn] 24Jul2015 


Dynamics at the Many-Body Localization Transition 


E. J. Torres-Herrera 1,2 and Lea F. Santos 1,3 
1 Department of Physics, Yeshiva University, New York, New York 10016, USA 
2 Instituto de Fi'sica, Universidad Autonoma de Puebla, Apt. Postal J-48, Puebla, Puebla, 72570, Mexico 
3 ITAMP, Harvard-Smithsonian Center for Astrophysics, Cambridge, Massachusetts 02138, USA 

(Dated: My 27, 2015) 

The isolated one-dimensional Heisenberg model with static random magnetic fields has become paradigmatic 
for the analysis of many-body localization. Here, we study the dynamics of this system initially prepared in a 
highly-excited nonstationary state. Our focus is on the probability for finding the initial state later in time, the 
so-called survival probability. Two distinct behaviors are identified before equilibration. At short times, the 
decay is very fast and equivalent to that of clean systems. It subsequently slows down and develops a powerlaw 
behavior with an exponent that coincides with the multifractal dimension of the eigenstates. 

PACS numbers: 72.15.Rn, 71.30.+h, 05.30.Rt, 75.10.Pq 


I. INTRODUCTION 

The metal-insulator transition has been at the forefront of 
physics research since Anderson’s seminal paper 01 . As a 
result of quantum interference, the wavefunctions of a disor¬ 
dered noninteracting system can become exponentially local¬ 
ized in configuration space. The phenomenon has been ex¬ 
perimentally observed in different setups, more recently with 
Bose-Einstein condensates 001. A proposal for an exper¬ 
iment with ultracold atoms in a two-dimensional geometry 
also exists 0]. Lattice models, such as the Anderson tight- 
binding and the powerlaw random banded matrix (PRBM) 
models 0H, have been extensively employed in the analy¬ 
sis of the Anderson metal-insulator transition. At criticality, it 
was found that the eigenstates exhibit multifractal features |5]. 

The inverse participation ratios, IPR“ = \C%\ 2q , con¬ 

tain information about the structure of the eigenstates | tf> a ) = 
]T n written in the basis vectors \<p n ) of the configura¬ 

tion space. In particular, IPRj measures the level of delocal¬ 
ization of the eigenstates in the chosen basis 0. At the An¬ 
derson transition, the probability amplitudes C“ display large 
fluctuations and IPR" shows anomalous multifractal scaling 
with respect to the system size I5ll8i- [l2ll . 

(IPR“) ~ U~ {q ~ 1)D \ (1) 

where (.) denotes the average over an ensemble of realizations 
and eigenstates, J\f is the dimension of the Hamiltonian ma¬ 
trix If 1311, and l) q represents the generalized dimension. Multi- 
fractality is reflected by the nonlinear dependence of the gen¬ 
eralized dimension on q. In contrast, D q = d in the metallic 
phase, where d is the system dimension, and D q = 0 in the 
insulating phase. Experimentally, multifractality has been ob¬ 
served in disordered conductors iHl and in systems with cold 
atoms El El- Recently, new studies have led to the con¬ 
clusion that multifractal correlations are not exclusive to the 
critical point of the Anderson-transition. In disordered sys¬ 
tems, they are present away from criticality El and even in 
extended states j 18]. They are also found in the ground states 
of clean systems |19fl. 

Studies of the dynamics of noninteracting systems at the 
metal-insulator transition have shown that the Loschmidt 


echo El . the survival probability Ell - El . and the spreading 
of wavepackets B23L l24fl at the mobility edge exhibit a pow¬ 
erlaw behavior, where the exponent coincides with the gen¬ 
eralized dimension for q = 2. The generalized dimension of 
the eigenstates, D 2 , is extracted from Eq. 0 by performing a 
scaling analysis of (IPR 2 ). However, in studies of dynamics 
the main interest is on the generalized dimension associated 
with the initial state \(f> no ) and denoted by D^. The latter is 
obtained from a scaling analysis of the level of delocalization 
of the initial state with respect to the energy eigenbasis, that is 
the analysis of 

(IPR 2 0 ) ~ Af ~° 2 , (2) 

where |(/>„„) = C“ 0 \tp a ) andIPR” 0 = J2 a \ C nX■ when 

investigating localization in real space, the initial state usually 
corresponds to a basis vector of the configuration space. We 
also note that the two generalized dimensions above have been 
shown [2l to be related through the expression D 2 = dD 2 - 

A natural question following this brief summary of the 
Anderson localization is what happens to the above findings 
when interaction is included. It had been conjectured al¬ 
ready in fl j . 1 2 61 and then confirmed with perturbative argu¬ 
ments E3l2Ml and rigorously [29|] that localization may per¬ 
sist. Studies about many-body localization (MBL) have re¬ 
cently boomed EoHEIl . The interest in the subject is in part 
motivated by the access to new experimental tools, such as 
cold atoms in optical lattices |53j, that can be used to corrob¬ 
orate theoretical predictions. Among the latter, we find works 
about the location of the critical point in disordered spin-1/2 
chains Ell - [34ll . analysis of the relation between the distribu¬ 
tion of the wavefunction coefficients and the onset of localiza¬ 
tion 1 32Sl . various efforts to identify the quasi-local integrals 
of motion in the MBL phase i47l450n . and descriptions of the 
evolution of the entanglement entropy 136, 37]]. few-body ob¬ 
servables GM], and the Loschmidt echo 1441. 

Our goal in this work is to characterize the evolution of iso¬ 
lated disordered systems with interaction from very short to 
very long times. Since MBL is a dynamical transition, identi¬ 
fying general features of the dynamics of interacting systems 
is essential for the further developments of the field. Moti¬ 
vated by the results for noninteracting systems, our focus is 
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on the decay of the survival probability and its relationship 
with the onset of multifractal states. 

We consider a one-dimensional (ID) disordered spin-1/2 
system and analyze the evolution at different time scales of 
both the survival probability and the time-averaged survival 
probability. At short times the decay is very fast and similar 
to that of clean systems. Afterwards, the decay slows down 
and shows an anomalous powerlaw behavior. The exponent of 
this algebraic decay coincides with D^. At very long-times, 
the decay eventually saturates to (IPRj 0 ). 


II. MODEL AND BASIS 


We investigate the ID isotropic Heisenberg spin-1/2 sys¬ 
tem with two-body nearest-neighbor interaction, L sites, and 
periodic boundary conditions. The Hamiltonian is 

L 

H = J2 [hkSl + J (stsi +1 + stst +l + S* k SZ +1 ) . 

fc=1 

(3) 

Above h = 1, S k ’ v,z are spin operators, and J = 1 sets 
the energy scale. Random static magnetic fields act on each 
site k, the amplitudes hk being random numbers from a uni¬ 
form distribution [—h,h\. The total spin in the z-direction, 
<S Z = ^2 k SI, is conserved. We work with the largest sub¬ 
space, S z = 0, of dimension M = L\/(L/ 2)! 2 . Localization 
in this symmetry sector guarantees localization in smaller sec¬ 
tors. 

The dependence on h of the level statistics and of the level 
of delocalization of the eigenstates of H ([3]) has been studied 
for at least a decade When h = 0, the system 

is analytically solvable with the Bethe ansatz. If all the triv¬ 
ial symmetries of the Hamiltonian are taken into account, one 
verifies that the level spacing distribution of neighboring lev¬ 
els is Poisson. In addition to the total spin in the z-direction, 
the other symmetries of the isotropic model at S z = 0 are: 
translational invariance, parity, spin reversal, and conserva¬ 
tion of total spin. 

As h increases from zero, the level spacing distribution 
eventually becomes Wigner-Dyson, indicating a transition to 
the chaotic regime. The value of h at which level repulsion 
becomes evident decreases as the system size increases. In 
parallel, the level of delocalization of the eigenstates in real 
space increases substantially. The presence of disorder breaks 
the additional symmetries mentioned above and if the disorder 
is weak the states can spread out significantly. 

As h further increases and becomes larger than the coupling 
strength, h > 1, the spreading of the eigenstates recedes and 
they become more localized in real space. The critical point 
for the transition to the MBL phase has been identified as h c = 
3.5 ± 1.0 in HI and h c = 2.7 ± 0.3 in l32[| . 

In any study of the structure of the eigenstates, the choice 
of basis is essential. Since here the goal is to investigate lo¬ 
calization in real space, that is the level of confinement of the 
spin excitations in the lattice, the natural basis is that of the 
configuration space, which we refer to as the site-basis and is 


also known as the computational basis. The site-basis vectors 
\(f> n ) correspond to states where the spin on each site either 
points up or down along the z-axis, such as | 


III. SURVIVAL PROBABILITY 


To study the dynamics of the disordered chain (0, we take 
as initial state a single site-basis vector, |\F(0)) = | cf> no ). This 
is equivalent to a quench, where the initial Hamiltonian is 
the Ising part of the Hamiltonian Yhk=i ^k ^k- i an d the fi¬ 
nal Hamiltonian is H 0. To quantify how fast the initial 
state changes in time, we concentrate on the behavior of the 
survival probability. 


Fit) 


(^(0)|e- lJ?t |T-(0)) 


Era 2 e-^‘ 
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(4) 


where E a are the eigenvalues of H and C“ 0 = {ip a \<fi no ) is 
the overlap of the initial state with the eigenstates \ip a ) of H. 
F(t) measures the probability for finding the system still in 
|\&(0)) at time t. 

The distribution in energy 

Pn 0 iE) = Y J \CZ 0 \ 2 KE-E a ) (5) 

a 

of the components |C“ | 2 of the initial state is often re¬ 
ferred to as local density of states (LDOS). If the envelope 
of this distribution is known, an analytical expression for 
F{t ) can be obtained from the Fourier transform, F{t ) ~ 

/ Pno (E)e-' Et dE. 

For strong quenches, that is when the initial and final 
Hamiltonians are very different, the envelope of p no (E) is a 
Gaussian with mean corresponding to the energy of the initial 
state, 

- E = (<^ol%no) (6) 


and width 

< = Era 2 ^ 2 -4 0 = E i<^i%n 0 >r. (?) 

a. n/no 

Notice that the width depends only on the off-diagonal ele¬ 
ments of the Hamiltonian matrix written in the site-basis and 
is therefore independent of the diagonal disorder. The Gaus¬ 
sian shape of the LDOS reflects the density of states, which 
for systems with two-body interaction is also Gaussian l55j]. 
In the absence of disorder, the envelope p no ( E ) is particu¬ 
larly well filled for initial states with energy e„ 0 near the cen¬ 
ter of the spectrum of // [56), 0|. Its Gaussian shape leads to 
the Gaussian decay F{t ) ~ exp(—a 2 0 f 2 ). This behavior may 
persist until saturation or be followed by an exponential (see 
Refs. [57j, Ht and references therein). 

In Fig. [I] we analyze the survival probability and the LDOS 
in the presence of disorder. The average of F(t) over differ¬ 
ent disorder realizations and different initial states is denoted 
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by ( F(t )). For each system size and each realization, we se¬ 
lect as initial states, only 10% of all the M site-basis vectors. 
They are the ones with energy e no closest to the middle of the 
spectrum of H. Since the density of states is Gaussian, the 
center of the spectrum contains the most delocalized states. 
Localization in this region assures localization in other parts 
of the spectrum. For each L, the total number of data points 
for the average, including initial states and realizations, adds 
up to 10 5 . 



FIG. 1: (Color online) Survival probability averaged over 10 5 data 
points for h = 0.5,1.0,1.5, 2.0, 2.7,4.0 from bottom to top (a) and 
LDOS for a single realization for the bottom panel h = 0.5 (d), the 
middle h = 1.5 (c), and the top h = 2.7 (b); L = 16. In (a): 
the dashed line indicates exp(—(<r^ 0 )f 2 ) for h = 0.5 and horizon¬ 
tal lines correspond to the saturation point, (IPRj 0 ). The envelopes 
(solid lines) of the distributions in panels (b), (c), and (d) are Gaus- 
sians with center e no [Eq. ©I and width a no [Eq. I©]. 


Figure [I] (a) displays ( F(t )) for different values of h. The 
initial decay is very fast until t ~ 2. For small disorder, the 
initial evolution is purely Gaussian. As h increases, the inter¬ 
val of the Gaussian decay shrinks until only the quadratic part 
persists, ( F(t )) ~ 1 — (cr„ 0 )f 2 . This is followed by a possible 
exponential behavior, the time interval being too short for cer¬ 
tainty. After the initial fast evolution, oscillations appear. The 
time interval of these oscillations as well as their amplitudes 
increase with the disorder strength. 

The oscillations eventually fade away and give place to a 
powerlaw decay with exponent < 1. The initial state finds 
new channels that give continuation to its evolution. The cou¬ 
plings at higher order in perturbation theory become gradually 
effective. 

The 1 o n g - ti mep o w e r law behavior reflects the onset of mul¬ 
tifractal states 1121 M2 311 . Our results indicate that multifractal 
many-body states can occur even at small h. As the disorder 
increases and the eigenstates become less extended, the pow¬ 
erlaw exponent naturally decreases. For L = 16, the decay 
after the oscillations is hardly noticeable for h > 4. 

At very long times, the decay eventually saturates. The 
saturation point is derived from the infinite time average of 
Eq.©, 


(F(t -» oo)) 



(IPR”°). 


The value of this infinite time average naturally increases with 
the disorder strength. 

Figures Q](b), (c), and (d) display representative LDOS for 
three values of h. The widths of the three distributions are 
equivalent, because according to Eq. ©. a no does not depend 
on the disorder strength. This explains the indistinguishable 
initial decay for all curves in Fig.Q](a). 

At small h [Fig. H|(d)], the Gaussian envelope of the distri¬ 
bution is still well filled, indicating a very delocalized initial 
state. This is independent of the realization, provided e no be 
near the center of the spectrum. As the disorder increases, 
the multifractal structures of the eigenstates spread to larger 
scales and the coefficients C“ o fluctuate strongly. As a result, 
the LDOS becomes more sparse [Figs. |T] (b) and (c)], justify¬ 
ing the oscillations and subsequent powerlaw decay in Fig. |T| 
(a). The oscillations are due to the small number of states en¬ 
ergetically accessible to the initial state in low order of pertur¬ 
bation theory; a number that decreases as h increases. These 
oscillations are not random fluctuations that can be averaged 
out with enough realizations, as those at very long times. They 
are connected with the approach to the MBL phase and the on¬ 
set of quasi-integrals of motion l l47l - l50ll . 

In Fig. [2] we analyze the survival probability for different 
system sizes and four values of h. The strengths of the dis¬ 
order are small in Figs. U (a), (b), while in Figs. [2] (c), (d), 
they coincide, within errors, with the critical point h c ob¬ 
tained in Refs. lf3lL [32ll . The fast evolution for t < 2 is sepa¬ 
rated from the later powerlaw decay either by a small plateau 
(a) or by visible oscillations (c), (d). The scope of the pow¬ 
erlaw behavior increases with system size and with disorder 
strength [compare the time where saturation takes place in (a) 
with the time in (c), for example]. This suggests that for very 
large L the algebraic decay may persist for h > h c . 



t t 


FIG. 2: (Color online) Survival probability averaged over 10 5 data 
point for h = 1.0 (a), h = 1.5 (b), h = 2.5 (c), and h = 2.7 (b) 
for L = 10,12, 14,16 from top to bottom. Dashed lines give t~° 2 , 
where D 2 = 0.99 (a), 0.84 (b), 0.36 (c), and 0.30 (d). 

The dashed lines in Fig. [2] correspond to an algebraic de- 
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cay described by the generalized dimension, ( F(t )) oc t ° 2 . 
As shown in Fig. [3] D 2 is extracted from the best linear fit to 
ln(IPR 2 °) vs In Af for L = 8,10,12,14,16. The error bars 
are standard deviations over 10 5 different values of IPR’',' 0 for 
each L. At small disorder, the error bars are small. As h ap¬ 
proaches h c , the dispersion of the values of IPR.; 0 and there¬ 
fore the uncertainty in the value of D 2 increases. At very large 
disorder, the error bars decrease again. 
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FIG. 3: (Color online) ln{IPR 2 °) vs In M (dark circle) and In IPR‘ 2 yp 
vs lnA/” (light square) for h = 1 (a), h = 1.5 (b), h = 2.0 (c), 
h = 2.7 (d), h = 3.2 (e), and h = 4.0 (f). Error bars are standard 
deviations over 10 5 values of IPR 2 0 (dark color) or of In IPR 2 0 (light 
color). 


For small disorder, h < 1, the system is still close to the 
metallic phase and the decay is diffusive, D 2 ~ 1. In this 
case, the exponent of the numerical powerlaw decay agrees 
extremely well with D 2 when the system size is large [see 
Fig. [21(a)] . As h increases, D 2 decreases, but not as fast as the 
numerical exponent. For h = 1.5 [Fig. [2] (b)], the agreement 
between the numerical curve and ( F(t )) oc t ~ D2 is not very 
good anymore. 

In the vicinity of the critical point. Figs. [2] (c) and (d), os¬ 
cillations are seen approximately in the same time interval of 
the algebraic decay of Fig. [2] (a). The generalized dimension 
is now Z ?2 < 1/2 and it agrees well with the rate of the damp¬ 
ing of those oscillations, while the powerlaw decay appears 
now latter in time. As L increases, the amplitudes of the os¬ 
cillations decrease and the slope of the subsequent powerlaw 
decay becomes more pronounced and closer to t ^° 2 . It is thus 
plausible to expect that for very large system sizes, D 2 might 
be able to capture the algebraic decay also for large disorder. 
This expectation is further supported by the results below for 
the time-averaged survival probability. 


IV. TIME-AVERAGED SURVIVAL PROBABILITY 

In the analysis of the dynamics of noninteracting systems at 
the mobility edge B20I423I1 . the commonly employed quantity 
is the time-averaged survival probability, which smoothes the 
fluctuations in ( F(t)). It is defined as, 

C(t) = - t J (F(t))cLt. (8) 

To reduce also the fluctuations in the values of IPR 2 °, one 
often deals with the so-called typical inverse participation ra¬ 
tio, I PR 1 ; 1 ’ = exp((lnIPR 2 °)). The scaling analysis of [PR' 2 yp 
gives Df, as shown in Fig. [3] The error bars for I PR 1 / 1 ’ in 
that figure are, of course, smaller than those for the regular 
IPR 2 0 , since now we deal with the dispersions in the values 
of In IPR 2 0 . 



FIG. 4: (Color online) Time-averaged survival probability for h = 
1.0 (a), h = 1.7 (b), h = 2.0 (c), and h = 2.7 (d) for L = 
10,12,14,16 from top to bottom. Dashed lines gives t 2 , where 
5 2 yp = 0.99 (a), 0.87 (b), 0.72 (c), and 0.42 (d). 

In Fig.0] we compare C{t ) with tr r> ^'. When the system 
is still close to the metallic phase, as in Fig. [4] (a), the decay 
of C(t ) is smooth all the way to saturation and in excellent 
agreement with t~ 2 , especially for L = 16. 

As the disorder increases, the powerlaw exponent de¬ 
creases, but the short-time dynamics does not change much. 
This creates an abrupt contrast between the two time scales, 
resulting in a visible elbow [see Fig. [4] (c) and (d)]. As h in¬ 
creases, we also notice that the time interval for the agreement 
between the algebraic decay of C(t) and A shortens and 
starts later in time (compare the four panels). Yet, for a fixed 
disorder, the agreement also improves with L. This indicates 
that for system sizes larger than available for exact diagonal- 
ization, D 1 ^ should be able to describe the powerlaw decay 
for long times, even when the disorder is strong. 
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V. POWERLAW EXPONENT AND SYSTEM SIZE 

The exponent of the powerlaw decay of F(t) contains im¬ 
portant information about the system: 

(i) Because it coincides with the generalized dimension D 2 , 
^indicates the level of delocalization of the initial state. Since 
D 2 ~ D 2 , as suggested in ||25j and confirmed by us for our 
model, the powerlaw exponent also manifests the level of de- 
localization of the eigenstates. 

(ii) It gives information about the correlations between the 
components |C“ | 2 , because the algebraic decay implies that 

lUeSlsI, 

<^C0> = (j2\ C no\ 2 \ C no\ 2 ^~ Ea)t ) 

\a,P / 

U}= E p—E a ^ f Kt -D2' (9) 

J ~ OO 

When the eigenstates, and consequently the initial state, are 
extended and thus similar to random vectors, as it happens 
in the chaotic domain (h < 1) for states close to the mid¬ 
dle of the spectrum, the components 16'" . | 2 are uncorrelated 
random numbers. In this case IPR 2 0 oc Af and the dynam¬ 
ics is diffusive ( D 2 ~ 1), as obtained also in (34)]. As the 
disorder increases, the states become multifractal; the com¬ 
ponents | C“ 0 1 2 show large fluctuations and become gradually 

more correlated, so IPR 2 0 oc Af° 2 with D 2 < 1, resulting in 
a subdiffusive dynamics. The limited spreading of the initial 
state quantified by D 2 reflects, as made explicit by Eq. 0, the 
level of correlations between the components |C“ 0 1 2 . 

Figure0(a) shows how D 2 and l) l J p depend on the disorder 
strength. In parallel with the standard deviations in Fig. [3} the 
error bars are larger for D > than for We show only the 
latter to simplify the figure. Within errors, the two generalized 
dimensions coincide. As h increases, the number of states that 
contribute to the evolution of the initial state shrinks and the 
generalized dimensions decrease. The decay is evident for 
1 < h < 4 and it becomes extremely slow afterwards. We 
avoid an analysis of what happens for h > 4, because for the 
very small system sizes L = 8,10 we actually see an approach 
to on-site localization. 

Even though it is not clear at this point how to identify the 
MBL critical point from Fig. 0(a), a comparison with previ¬ 
ous studies is instructive. The values of h for the mid-point be¬ 
tween a metal and an insulator, that is D^ ~ 1/2 (h ~ 2.5), 
and for the inflection point of the fitting curve (h ~ 2) are not 
too far from the critical points found in 011 13211 . In addition, 
the point of an almost halt in the decay of the values of the 
generalized dimensions, h ~ 4, is very close to the critical 
point h c ~ 3.7 obtained in Ref. m for states that, as in our 
case, live close to the middle of the spectrum. 

We expect D 2 and to get closer for scaling analysis 
performed with larger system sizes than the very few ones 
now available for exact diagonalization. As L increases, they 
should also better agree with the powerlaw exponent of F(t). 
These claims find support already in the results for L = 14 



FIG. 5: (Color online) D 2 (dark circle) and D'l v (light square) vs 
disorder strength (a) and the distribution of ln(IPR() 0 /(i) for h = 2.7 
and L = 8, 10,12, 14,16 from highest to lowest peak (b). In (a): 
solid lines are fitting curves. Only the error bars for D 1 ^ are shown. 
They are smaller than those for D 2 , in accordance with the standard 
deviations in Fig. [3] 


and 16. If the powerlaw behavior is indeed to be described 
by the generalized dimension, then the algebraic decays for 
different system sizes must coincide. In Figs. U and [4] the 
slopes are visibly different for small L’s, but they get closer 
for L = 14 and 16. This suggests that the scaling analysis 
should become more accurate for sizes L > 14. 

Figure 0(b) endorses the proximity of the results for L = 
14 and 16. It shows the distribution of the inverse partici¬ 
pation ratios. IPRIJ 0 fluctuates with disorder realization and 
initial state. However, the validity of Eq. (0 presupposes that 
D 2 does not depend strongly on what is used on the left side 
of that equation, whether it is (IPR?? 0 ), IPR 1 / 11 , or the most 
probable value of IPRj 0 . This implies that the distribution 
of IPRj 0 normalized to its median /< must have a scale in¬ 
variant shape mu. As seen in Fig. 0(b), the distribution 
P[ln(IPR 2 °//r)] broadens considerably from L = 8 to 12, but 
the shapes are similar for L = 14 and 16. 

In noninteracting disordered systems described by the pow¬ 
erlaw random banded matrix, numerical evidence for the scale 
invariance of PpiiIPR^ 0 ] was achieved Hill already for Af > 
300, in contrast with the Af > 3000 needed here. The exis¬ 
tence of more correlations between the matrix elements of our 
system when compared to random matrices may justify such 
large difference. The number of nonrandom elements in the 
Hamiltonian matrix of Eq. 0 is much larger than in noninter¬ 
acting systems, such as those described by the tight-binding 
model or the powerlaw random banded matrix. 


VI. CONCLUSION 

We studied the dynamics of an isolated disordered ID 
Heisenberg model as it approaches the MBF phase. The anal¬ 
ysis was based on the entire evolution of the survival probabil¬ 
ity F(t ), from t = 0 to t —> 00 , for initial states corresponding 
to site-basis vectors. F(t) is one of the simplest quantities that 
can reveal the multifractality of the eigenstates. It also appears 
explicitly in the evolution of observables |[57l l. 

The dynamics of clean and disordered interacting systems 
is comparable at short times. For both, the Gaussian decay 
rate of F(t) coincides with the width of the FDOS. In the 
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presence of disorder, the LDOS gets sparse, reflecting the re¬ 
duced number of states participating in the dynamics and the 
multifractality of the eigenstates. As a result, the behavior of 
F(t) at long times becomes powerlaw. 

The exponent of the powerlaw decay coincides with the 
generalized dimension Di. This finding establishes a paral¬ 
lel with previous works about the dynamics of noninteracting 
systems at criticality and may help advance our understanding 
of transport properties in interacting systems. It also implies 
that from F(t), one can infer the level of delocalization of the 
initial states and eigenstates, as well as the correlations of their 
components. This is advantageous, since numerical methods 
other than exact diagonalization are available for studying dy¬ 


namics, which gives access to larger system sizes. The dy¬ 
namics can also be studied experimentally with quantum sim¬ 
ulators. 
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